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Twisted-mass lattice QCD with mass non-degenerate quarks * 
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The maximally twisted lattice QCD action of an SUf{2) doublet of mass degenerate Wilson quarks gives 
rise to a real positive fermion determinant and it is invariant under the product of standard parity times the 
change of sign of the coefficient of the Wilson term. The existence of this spurionic symmetry implies that 0(a) 
improvement is either automatic or achieved through simple linear combinations of quantities taken with opposite 
external three-momenta. We show that in the case of maximal twist all these nice results can be extended to the 
more interesting case of a mass non-degenerate quark pair. 
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1. Introduction 

It has been shown in ref. that in lattice QCD 
(LQCD) with mass degenerate SUf{2) doublets 
of Wilson fermions it is possible to improve the 
approach to the continuum limit of correlation 
functions of gauge invariant multiplicative renor- 
malizable (m.r.) operators by taking the arith- 
metic average (Wilson average - WA) of pairs of 
correlators computed in theories regularized with 
Wilson terms of opposite sign and identical values 
of the subtracted (unrenormalized) lattice quark 
mass, — Mq — Ma- with A/q the bare quark 
mass. Equivalently one can take appropriate lin- 
ear combinations of the correlators computed in 
a given regularization (fixed sign of the Wilson 
term) but with opposite values of (mass av- 
erage - MA) . Improved hadronic masses and ma- 
trix elements can be similarly obtained by taking 
VFA's of the corresponding quantities separately 
computed within the two regularizations. 

To avoid the difhculties related to the nature 
of the spectrum of the Wilson-Dirac (WD) op- 
erator (here we are referring to the well-known 
problem of the existence of "exceptional config- 
urations" j2l3j ) twisted-mass lattice QCD (tm- 
LQCD) 15 should be used for the actual com- 
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putation of the correlators taking part to the av- 
erages. The fermionic determinant of tm-LQCD 
is, in fact, real and positive on arbitrary gauge 
backgrounds, as long as the quark mass is non- 
vanishing. All the nice cancellations of 0(a) 
terms that one finds in the standard Wilson case 
extend to tm-LQCD with mass degenerate quark 
doublets. 

Peculiar simplifications occur if the twisting an- 
gle is taken to be equal to 7r/2 (maximal twist). 
This choice is particularly useful for applications 
because essentially all interesting physical quan- 
tities (e.g. hadronic masses and matrix elements 
of m.r. operators between states with vanishing 
three- momenta) can be extracted from 0(a) im- 
proved lattice data without making recourse to 
any WA. 

The main goal of this note is to show that in the 
maximally twisted case positivity of the determi- 
nant and 0(a) improvement without WA can be 
extended to encompass the more interesting sit- 
uation in which mass non-degenerate quark pairs 
are considered. This is a preliminary, necessary 
step if one wants to set up a realistic computa- 
tional scheme for operator matrix elements. 

For future applications it is important to stress 
that the approach we propose is very flexible as 
it allows to regularize different flavours with dif- 
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ferently twisted Wilson terms. As discussed in 
detail in ^ , it is possible to exploit this freedom 
in order to improve the chiral behaviour of lat- 
tice correlators to the point of (hopefully) killing 
all the unwanted "wrong chirality mixings" that 
affect in the standard Wilson |H] and twisted- 
mass |4I7| regularizations of QCD the construc- 
tion of m.r. operators. Of special importance is 
the application of this new strategy to the con- 
struction of the m.r. operators which represent 
the CP-conserving, AS" — 2 and AS" — 1 effective 
weak Hamiltonian on the lattice. 

The presentation of the material is divided in 
two parts. In the first part (sect.[5J) we spell out 
in detail the formulae relevant for mass degen- 
erate quark pairs with the choice ±7r/2 (maxi- 
mal twist) for the twisting angle. In the second 
part (sect. we discuss how one can deal with 
mass non-degenerate pairs of quarks without loos- 
ing 0(a) improvement, or the positivity of the 
fermion determinant. The latter property, which 
is obvious in the case of mass degenerate dou- 
blets, is crucial for being able to actually carry 
out numerical Monte Carlo simulations. Conclu- 
sions can be found in sect. ^ In Appendix A 
we prove the renormalizability of the maximally 
twisted fermion action both in the case of mass 
degenerate and mass non-degenerate quarks. The 
proof of the positivity of the fermion determinant 
in the non-trivial case of mass non-degenerate 
quarks is given in Appendix B. 

2. A summary of tm-LQCD @ ±tt/2 

The twisted- mass lattice action of an SUf{2) 
flavour doublet of mass degenerate quark has the 
form 



5^^[V',^,C/]=a4^^(x) 

X 

+e-'"^="^W/cr(r)+m,]v(a;) 
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(1) 



where 



(2) 
(3) 



and M„{r) = —Mcr{—r) is the critical quark 
mass. We wrote the action in what is usually 
called the "physical basis" jT], where m^, is taken 
to be real (and positive) ^ . By undoing the twist- 
ing of the Wilson term and bringing it fully to 
the mass term, it was proved in ref. ^ that the 
parameter Mcr(r) appearing in the action is 
equal to that of the corresponding standard Wil- 
son theory. 

In this note we will concentrate on the maximal 
twist case oj = tt/2, where eq. becomes 



X 
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(4) 



The choices uj = ±7r/2 are particularly useful be- 
cause it can be proved ^ that, despite the fact 
that the theory is not fully 0(a) improved, can- 
cellation of 0(a) effects in quantities of physical 
interest (like energies and operator matrix ele- 
ments) is either automatic or obtainable with no 
need of any WA. The proof of this statement 
is sketched in sect. 12.21 0(a) ambiguities in the 
knowledge of Mci{r) do not spoil any of the above 
results. 

Furthermore it is important to observe that at 
bj — ±7r/2 the critical WD operator (i.e. the oper- 



ator in square parenthesis in eq. I@J with : 



0) 



is anti-Hermitian, so its spectrum is purely imag- 
inary. This means that the full WD operator 
cannot have any vanishing eigenvalue as soon as 
Tfiq ^ 0. This is also evident from the explicit 
expression of its determinant, Pf.d, which takes 
the remarkably simple form 



Pf.d = det Ql^Q. 



(5) 



where Qcr = Q\i is defined by eqs. (0) and H69|l . 

In Appendix A we prove that the fermionic ac- 
tion is "stable" under radiative corrections, 
in the sense that the unbroken (possibly spu- 
rionic) symmetries of the action prevent radia- 
tive corrections from generating extra indepen- 
dent operators (of dimension < 4), not already 

^Unless differently stated, wc employ here the notations 
of ref. m. 

^Whatever we say for ui = 7r/2 also holds for ui = — 7r/2. 
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present in which one would need to include 
for renormalizability. Once Mcr{r) has been set 
to the appropriate value, the continuum limit is 
approached as usual by rescaling and nig ac- 
cording to the chosen renormalization conditions. 

2.1. Non-singlet Ward-Takahashi identi- 
ties 

In this section we collect the expressions of the 
renormalized current and quark density operators 
entering the flavour non-singlet Ward-Takahashi 
identities (WTI's) associated with the action Q). 

Renormalized vector and axial currents can be 
taken to be 



v.} 
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(6) 



where to make contact with known quantities the 
(finite) renormalization constants, Zy and Za, 
introduced above are those for the local vector 
and axial currents of standard Wilson fermions, 
respectively. Notice the switch between Zy and 
Za for the currents with flavour h — 1,2, due 
to the presence of the factor 75X3 in front of the 
Wilson term in eq. ^ (see also the comment at 
the end of this section). 

With reference to eqs. ©, the non-singlet 
WTI's with the insertion of the renormalized 
(multi-local) operator 0{y) {y = {yi,i = 
1, . . . , 71} x) take the expected form {b = 1, 2, 3) 



{dlV'^{x)d{y)) 



{ dlAiix) - 2m,P\x) 0{y)) 



(r,m,) 



0(a) 



(r,i 



with 



rhq = Zp^ruq , 



(7) 

= 0(a) (8) 
(9) 



provided one deflnes, in terms of bare quantities, 
the renormalized pseudo-scalar operators, P^, to 
be 



1,2 



(10) 



= Z.0 



V'y75V' + a ^ipp{amg)l . (11) 



In eq. Hll|) the dimcnsionless real coefficient, 
pp(amq), in front of the power divergent term 



admits a polynomial expansion in airiq. Similarly 
to Zy and Za, Zp and Zgo are the (logarithmi- 
cally divergent) renormalization constants of the 
non-singlet pseudo-scalar and singlet scalar quark 
density of the standard Wilson regularization, re- 
spectively. 

For completeness we also give the expression 
of the renormalized singlet scalar quark density 
operator, S^. By standard symmetry and dimen- 
sionality arguments, one finds 



S" = Z, 



+ a ^mq/dgo (arrig)! 



(12) 



where ^50 is a dimcnsionless real coefficient with 
an even polynomial dependence on amq. This 
parity property follows from the invariance of 
the action under the spurionic symmetry 

^^1/2 ^ i^g ~^ -'^q) (eq- IHI))- Formula lO is 
rather interesting because it shows that the chi- 
ral order parameter is only affected by an m^/a^ 
power divergence, analogously to what happens 
with Ginsparg- Wilson fermions. 

The equations of this section can be obtained 
by specializing to the case = w = 7r/2 the 
corresponding formulae of ref. A direct and 
practical method to deduce them is to proceed 
in the following way. One starts by rotating the 
fermionic fields in the action Q so as to have the 
Wilson term in the standard form (i.e. with no 
factor 75T3 in front of it). At this point one can 
straightforwardly adapt to this instance the anal- 
ysis of the structure of the mixing pattern of the 
chiral rotation of the Wilson term spelled out in 
detail in ref. H] (see also ref. fS!). In a mass inde- 
pendent renormalization scheme renormalization 
constants can be computed in the massless limit 
{niq — 0). For this reason it seemed to us a natu- 
ral choice to keep for them the names they would 
have in the standard Wilson regularization. The 
final step of this method requires to rotate all 
the quark fields back into the physical basis we 
started from. 

2.2. 0(a) improvement 

It was proved in ref. that the invariance of 
the action Q under 

nf = 7^5 X (r ^ -r) X (m, ^ -m,) , (13) 
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(14) 



7^5 X , 



with 72.5 and Vd defined in eqs. and ||0i^ 
respectively, implies the vahdity of the formula 



(u) WA 



(O) 



(O) 



cont 
(m<j) 



(15) 



where O is any gauge invariant m.r. (multi-local) 
operator. 

From H15|l one can prove a number of phys- 
ically interesting improvement formulae. To see 
this let us introduce the eigenstates \h, n, k)||^|^ ^ 
{h and k represent the set of quantum numbers 
and three-momenta characterizing the state and 
n the excitation level) of the transfer matrix, 
T{uj,r,mq), with eigenvalues Eh^n(}^',^,'i',mq) ^. 
The eigenvalue equation reads 



T{u),r, mq)\h,n, k) 



(UJ) 



(r,m,) 



(16) 



In the notations of ref. ^ one gets the formulae 
£^h,„(k; w,r,TOq) -f (k; tj, -r, m,) = 



2i?n*(k;™,) + 0(a2)^ 



(/i,n,k|B|/i',?i',k') 



(17) 



+ (r ^ -r) (18) 



= 2(1 (w, r)(/i, n, k\B\h',n', k'} -t- 0(0') , 

(m,) 

where B is a gauge invariant m.r. local opera- 
tor. It is important to remark that in the whole 
argument about 0(a) improvement of VFA's the 
twisting angle w is a totally inert label. 

The interesting observation is that at lu — Tr/2 
the second term of the VFA's in eqs. H17(l and H18|l. 
i.e. the quantities evaluated with Wilson parame- 
ter — r, can be rewritten in terms of closely related 

* Actually the transfer matrix has been constructed only 
for \r\ = 1. For \r\ < 1 one should make reference to 
its square For simplicity we do not do that in the 

following. 



quantities but evaluated with the value r of the 
Wilson parameter. As a consequence, 0(a) im- 
proved estimates of energies and matrix elements 
can be obtained without having to average results 
from simulations with lattice actions differing by 
the sign of the Wilson term and critical mass. 

Before coming to the proof of this statement, 
we recall from ref. ,lj, that at generic values of u) 
the product 

Vx{uj^-uj), (19) 

where V is the physical parity operator (xp = 
(-x,t)) 




Uo{xp) , 
70-0 (a;p) 



ak),k^ 1,2,3 



(20) 



is a spurionic symmetry of the tm-LQCD ac- 
tion ((TJ. Since [V x {ui ^ = 1, (multi- 
local) m.r. operators can always be taken to have 
a definite parity, which can be read off from the 
formula 



{o^'Hix.})) 



(r,mg) 



^ {-If {0^^\{x,p})) 

This relation entails the possibility of defining a 
notion of parity for the eigenstates of the transfer 
matrix 1 . One can, in fact, prove the validity of 
the following equations 

VT{uj, r, mq)V = T{—lo, r, m, 
T]h.n\h,n, 



V\h,n,k) 

(r,m,) 



£^/i,n(— k; 



-uj, r, niq) . 



(21) 

(22) 
(23) 



where V is the representative of the parity oper- 
ation on the Hilbert space of states of the theory 
and r]fi,n (Vh n ~ ^) what we will call the par- 
ity of the state \h,n,'k)\''^^^ y Indeed rjh.n can 
be taken to be an w-independent integer coincid- 
ing with the physical parity label of the corre- 
sponding continuum state. For details we refer 
the reader to Appendix F of ref. 

It is also immediate to recognize that the ac- 
tion ^ goes into itself under the transformation 
(we recall that uj is defined mod 27r) 

(r —r) X {uj ^ UJ ±tt) . (24) 
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If in particular we set, say, lo = 7r/2, the tm- 
LQCD action takes the form and the invari- 
ance (|24|) becomes 



(25) 



Either by inspection or by observing that, owing 
to 1211), ifT^ is equivalent to 



(26) 



we conclude that the action I^J is invariant un- 
der (|26|l . This means that at uj — tt/2 eq. (|16|l 
and equations from (|21|l to H23I) can be rewritten 
in the form (for short we drop the twisting angle 
label) 



T{r, mq)-e 



\h,n,k) 



(r,m,) 



VT{r, mq)V = T{-r, ruq) , 



V\h,n,'k) 



(r\m,) 



£^ft,n(k;r,TOq) = £;,i,„(-k; -r, m,) . 



= 0(27) 

(28) 
(29) 
(30) 



Using the last two equations it is now possible 
to cast at w = 7r/2 the WA's (0 and (HTjl in 
a form in which all the relevant lattice data are 
extracted from simulations carried out with the 
action I^J and a given fixed value, r, of the Wilson 
parameter. One gets, in fact, the formulae 



(/i,n,k|B|/i',n',k') 



(r\m,) 



+vLh'n' {h, n, -'k\B\h', n , -k') 



(r,m,) 



= 2C|(r)(/i,n,k|B|;i',,^',k') +0(a2), (31) 

(m,) 

where ilhn,h'n' = Vhni-^y^'Vh'n' = ±1, and 

Eh,nO<-; r, niq) + Eh,n[~'^] r, niq) = 

= 2Ef,°^\k;mq) + 0{a'). (32) 

We observe that 77^^ is the product of the par- 
ities of the states |ft.,n, k) and |/i',n',k') (which, 
if not known, can be determined numerically as 
explained in ref. PP) times the parity, (— 1)^^, of 
the local m.r. operator B. 



Notice that in case the lattice matrix element 
(/i, n, k|_B|/i', n', k')|(r ,„ ) is invariant under in- 
version of all the external three-momenta (like 
when all three- momenta vanish) , the formula H31|l 
gets particularly simple. In fact, if 77^^ = 1, 
the lattice matrix element turns out to be auto- 
matically 0(a) improved, while if ry^^ = — 1 
eq. (|31|l implies that the 0(a) improved estimate 
of (/i, n, k|i?|/i', rt', k')|(^ „j ) is zero. This last re- 
sult is in agreement with what one expects in the 
continuum limit from parity invariance. 

An important instance of a quantity which is 
automatically 0(a) improved is -FV as it can 
be extracted from two-point correlators evaluated 
at zero external three-momentum. 

Eq. (|32|l taken at k = tells us that masses are 
automatically 0(a) improved. 

We end by observing that, although hadronic 
masses are automatically 0(a) improved at w = 
±7r/2, it is not possible to extract Mcr from nu- 
merical simulations with a discretization error 
which is better than 0(a) (unless the theory is 
fully 0(a) improved a la Symanzik). The rea- 
son is precisely that at cj = ±7r/2 the pion mass 
is automatically 0(a) improved, even if the crit- 
ical mass is known with an 0(a) error (see Ap- 
pendix D of ref. ^), thus it cannot be sensitive 
to 0(a) discretization uncertainties in the critical 
mass. With an analogous argument a similar con- 
clusion is reached if WTI's are used to determine 

3. Mass non- degenerate quarks 

The maximally twisted LQCD fermionic action 
of an SUf{2) pair of mass non-degenerate quark 
can be conveniently written in the form 



^?ii?][V',Vi,t/]=a4^Vi(a;)[7.V 

X 

-175X1 Wcr(r-) +mq + T3eq'^ip{x) , 



(33) 



where to keep the mass term real and flavour di- 
agonal we have used the matrix T3 to split the 
masses of the members of the doublet. Conse- 
quently the Wilson term was twisted with the 
flavour matrix ti. Our notations are such that 
TTiq and Sq are both positive. 
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What we said about the renormahzabihty of 
the fermion action in the mass degenerate case 
is valid also here and it is exphcitly proved in 
Appendix A. 

3.1. Non-singlet Ward-Takahashi identi- 
ties 

For mass non-degenerate quarks non-singlet 
WTFs take the continuum-like form {x ^ y) 



d;V^{x) + 2zi,S\x)\d{y)) 



(r,m,,eg 



{dlV^{x)0{y)) 



0(a) 



( dlAl{x)-2rn,P\x)\d{y)) 
dlAl{x)-2fn,P\x)\d{y)) 



(r,m,,c,) 



= 0(a) 
= 0(a) 

0(a) 
0(a) 



( dlAl(x)~2m,P^{x) + 
i,P\x)\d{y)) 



0(a). 



(r,m,,e,) 

if we make use of the definitions 

= Zv ^^7/. K = Za V^7a.75 f V' 

V^ = ZAi^l^,^^ Al = Zv 4^7^15 (34) 
= Za ^^7^ f ^ = ^^7^75^^' 



as the renormalized masses of the quarks in the 
doublet. 

For completeness we record the formulae 



5" = Zp 
- Zs 



ip^p + a ^mqPso{amq,aeq)li (43) 
?/'yV + a"^e9Ps(am,,ae,)a ,(44) 



where pso{amq, aeq) and ps{amq,aeq), as well as 
pp{amq,a€q), are dimensionless real coefficients 
that admit a polynomial expansion both in arriq 
and aCq. The mixing coefficients pp{amq) and 
Pso{amq), appearing in eqs. and (|f 2() in the 
mass degenerate case, coincide with the value at 
Eg = of the coefficients introduced in eq. (|35(l 
and (|43|l . respectively. 

3.2. 0(a) improvement 

The method for 0(a) improvement is just as in 
the mass degenerate case. The physical explana- 
tion of this fact is obvious: all 0(a) discretization 
effects come from the Wilson term which is of the 
same form in the actions Q and (up to a 
trivial flavour rotation). Formally, it is enough 
to observe that the action (|33|l is invariant under 
the transformations TZ^ x and 



7^' 



spND 



(45) 

7^5 X (r -r) X (m, -"?,,) x (e, -eq) , 



pb 



Zao 



V'y75V' + a ipp{amq,aeq)ll (35) 



ZP 75'/' : 



& = 2,3 



P" = Zs V'75V' 



(36) 
(37) 

(38) 
(39) 

(40) 

As expected, the above formulae turn into eqs. © 
to if we set = and perform the cyclic 
permutation of flavour indices 3 ^ 2 — > 1 ^ 3. 
The above WTFs allow us to identify 



Zpa V-y V' 



Zp^niq 



Zg ^€q . 



,(+) 



(q — Zp^niq 



iq — Zp^rriq — Zg ^tq 



(41) 
(42) 



as well as under the spurionic parity opera- 
tion H2t)|l . Invariance under TZ^^^^ and TZ^ x Vd 
allows to prove (by arguments analogous to those 
presented in ref. ^) that WA's of energies and 
matrix elements are free from 0(a) cutoff effects. 
In formulae (dropping the label lo — tt 12) 



Eh.n{^; r, niq, e,) -f £^h,„(k; -r, m,, e,) 



= 2i?;j°f (k;m„e,) + 0(a2). 



(46) 



(/i,n,k|B|/i',n',k') 



-I- (r — > -r) 



2C^(r)(/i,n,k|B|/i',n',k') 



0(a').(47) 



Using the symmetry of the action under 
the transformation (|26|l . one can then prove, 
along the lines of sect. 12.21 that averages of ener- 
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gies evaluated with opposite values of the three- 
momenta as well as the appropriate linear com- 
binations of matrix elements of m.r. operators 
taken between pairs of states with opposite ex- 
ternal three- momenta are not affected by 0(a) 
discretization errors. In formulae 



i?/i,„(k; r, mq,eq) + i?/i,„(-k; r, m,, e,) 



2£;£°f(k;™„6,)-|-0(a2) 



(48) 



{h,n,k\B\h',n\k') 
+Vhnh'n'{h,n,-k\B\h',n',-k') 
= 2Ci{r){h,n,k\B\h\n',k') 



+ 0{a^) .(49) 



3.3. The fermion determinant 

The fermion determinant associated with the 
action H33I) is real and strictly positive, provided 



ml > t\ 



(50) 



The proof, though elementary, requires some al- 
gebra and is presented in detail in Appendix B. 

It is important to realize that H50|) is not a triv- 
ial condition, because in terms of the renormal- 
ized quark masses, it implies the inequality 



,(-) 



(51) 



The latter leads to a rather stringent constraint 
on the (finite 
gets, in fact 



on the (finite) ratio ZpjZg^ if rriq^ ^ rriq \ One 



Zs 



> 1 



,(+) 



O 



(52) 



Since numerically one finds [lOlllj for the 

quenched ratio Zp/Zs a number somewhat 
smaller than 1, one might be worried that the 
inequality (|51|l is not fulfilled. On this issue a 
leverage can be offered by the choice of the value 
of r which can be any real number satisfying 
< |r| < 1. In perturbation theory one finds 
for the standard Wilson action to one-loop 
something like Zp/Zs = 1 ^ gif^Iif) + • . •, where 



I{r) — I{~r) > for < |r| < 1 and has a finite 
limit at r = 0. This suggests that decreasing |r| 
may increase Zp jZg . 

4. Conclusions 

In this note we have shown that tm-LQCD at 
io = ±tt/2 yields a particularly useful lattice regu- 
larization of the Wilson type for the gauge theory 
of mass non-degenerate fermion pairs. Following 
ref. we have proved that with this action es- 
sentially all physically relevant quantities can be 
evaluated with no 0(a) cutoff effects, while at 
the same time having the fermion determinant 
real and strictly positive, provided < < m^. 
Monte Carlo simulations are hence safely feasible. 
The difficulties of the standard HMC in the pres- 
ence of a flavour non-diagonal structure of the 
WD operator, I>nd (see eq. (IHTH . can be over- 
come e.g. by stochastically computing its deter- 
minant using algorithms of the multi-boson 
or PHMC 13, type. Such algorithms should now 
be based on some polynomial approximation of 

1/ dJjqIJnd and naturally allow a simple cor- 
rection of the employed polynomial approxima- 
tion. 

The approach we have discussed can be 
adapted |S] to the calculation of matrix elements 
of the CP-conserving, AS" = 2 and AS — 1 ef- 
fective weak Hamiltonian, in a way which we ex- 
pect will solve the problem of "wrong chirality 
mixings" in the construction of the correspond- 
ing m.r. lattice operators. 
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Appendix A: Symmetries of the action (|33|) 
and absence of counter-terms 

We want to show that the (spurionic) symme- 
tries enjoyed by the fermionic action H33|) ensure 
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that no extra operators (of dimension < 4) can be 
generated by radiative corrections, besides those 
already present. The discussion that foUows ap- 
phes in particular to the case = 0, implying 
the "stability" under radiative corrections of the 
action Q ^. 

The list of transformations which leave the ac- 
tion H33|l invariant include 

1) lattice gauge transformations, space-time 
translations and hyper-cubic rotations, as well 
as the vector transformation associated to 
baryon number conservation (none of these will 
play any special role in the discussion) 

2) three continuous (one vector, I^{9i), and two 
axial, P{92), I^{9^)) non-singlet spurionic trans- 
formations 



X^(6'i) X (m, -f e,r3 



i^i/2r 



e<;T3 e 



(53) 
(54) 



I\e:i) X (m, +e,r3 ^ 

^ g-*e375r2/2 ^^^^ ^ g^^g]g-*e375r2/2 ) ^ (55) 

where (notice the simplification occuring iih = 1) 
3) charge conjugation, C 



C : { 2p{x) "> i7o72V''^(a;) 

?A(a;) -jV'^(a;)7o72 



(56) 



4) the (anti-linear) reflection operation, 8. It 
can be defined either with respect to the time 
slice xq = (site reflection, Qg) or xq = a/2 
(link reflection, Qi). Its action on monomials of 
fermionic fields is 

e[/(C/)^(xi)...^(a:„)] = 

= r(e[[/])e[^(x„)]...e[V'(xi)], (57) 

where f{U) is a functional of link variables and 

Os/i[Ukix)] ^ u t{e,„x) , 

^The action 1331 at = is related to the action Jl} 
by the harmless flavour rotation if) — > exp (— iT27r/4)i/', 
i/i — > ?/iexp (ir27r/4). 



e./£[J7o(x)] = U^{e,„x - aO) , 



(58) 
(59) 



with 



9e{-xi,t) = {■x,~t + a) , 6l,(x,t) = (x, -t) (60) 

5) the pseudo-parity transformations 
'Pl/2 X i^h ^ -mq) , T'l X (e, ^ -e,) , P| , (61) 
where (xp = (— x, xq)) 



P. 



(62) 



C/o(a;) ^ ?7o(a;p) , 
Uk{x) ^ (72(xp -afc) ,fc= 1,2,3 
''/^ ' ) '(/;(x) ^ i75Ti7oV'(a;p) 
V'(x) ^ iV'(a;p)7o75n 

;7o(.t) ^C/o(a;p), 
^2,3.) J7fe(a;) f/2(xp - aft) ,A: = 1,2,3 .g^x 
■ ' ^/j{x) -> ir2,37oV'(a;p) 
V'(a;) -i^j{xphoT2,3 

6) the transformation TZ^ x P^, where 



^{x) -iIj{x)j5 

Uf,{x) Ul{-x - afi) 

'Dd : { i^{x) e3»'^/V(-a;) 

i^{x) e3"/2v;(-2;) 



(64) 



(65) 



7) the spurionic transformations (|26|l . (|45|l and 

Te x (r — > — r), where Te is the (euclidean) time- 
axis inversion. The action of Te on the elemen- 
tary fields of the theory is {xt = (x, —xq)) 



Uo{x) UlixT - aO) , 
Uk{x)^Uk{xT) , ft =1,2,3 
ip{x) -ypiQil^ixr) 
ip{x) ^{xT)lolb 



(66) 



In Tables 1 and 2 we list all the independent op- 
erators of dimension not larger than 4 that can- 
not appear in the action H33() . They are grouped 
in columns under the name of the corresponding 
"killing" symmetry. The reality properties of the 
various coefficients in the action l|33|l are fixed by 
the invariancc under the anti-linear reflection op- 
eration 6f (or Qs). 
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Table 1 

We list under the name of the relevant symme- 
try that forbids them all the operators of dimen- 
sion d = 4 and d = 3 that cannot appear in the 
density action (|^ . By FF we mean any lattice 
discretization of e^yApTr[F^yFAp]. 



Dim 


C 


d = 4 


^T757 Wi/j (T = 1, Ti, Ts) 


d 3 





to note the self-adjointness relation 
75T3£'nd75T-3 = D'Ib ■ (68) 

The proof of positivity is somewhat more in- 
volved. To proceed it is convenient to introduce 
the auxiliary self-adjoint operator 



75 



7 • V -I- Wcr{r) 



(69) 



in terms of which 2?f,nd can be written in the 
form 



2?F,ND = Det[(3c 



(70) 



In flavour space the operator in the r.h.s. of 
eq. H7()|) is represented by the 2x2 matrix 



ND e 



Table 2 

Continuation of Table 1 



Dim 




P'f X (e? ^ 


d = 4 


^pTj^tp (T = Tl,T2) 

FF 




d = 3 







The conclusion of this discussion is that the 
form H33(l of the action is preserved by radiative 



corrections and the bare parameters. 



m, 



and 



Eq, need only a purely multiplicative renormaliza- 
tion 
and e 



More details on the rcnormalization of niq 
q can be found in the text in sect. 13. ll 



Appendix B: Positivity of the fermion de- 
terminant 

In this Appendix we show that under the con- 
dition H50() the determinant of the WD operator 
associated with the action (|33|l . namely 



D 



ND 



J ■ V - ij5TiWcr{r) + niq + Tseq , (67) 



is real and (strictly) positive. To prove that 
2?F,ND = Dct[I?ND] is a real number it is enough 



Qcr + 75e9 



(71) 



•■q Qci ~ l5^q 

For the determinant of this operator one finds 
2?F,ND = Det[i:)ND] = 

= det[((3cr + 75eg)(Qcr ~ l5<^q) + m^q] = 

^det[Qi + ml-€l + 2eq^-V] = 
= det[Q2^, + TO^ - e^] ■ dct[l + l^qB] , (72) 
where we have introduced the definition 
i? = A7.VA, A = (Q2^ + m^-e2)"V2.(73) 

Since det[Q^j.-|-mg — e^] > 0, we only have to prove 
that the second factor in the last equality of H72|l 
is a positive quantity. To this end we observe that 
both 2?F,ND (from eq. (|7T|l ) and det[(5^, 
(by inspection) are even under ^ — e 
then conclude that 



We 



AF,ND(e9) = dct[l + l^qB] = AF,ND(-e9) . (74) 

This property is indeed enough to prove the pos- 
itivity of Af,nd and hence our thesis. Let us 
consider, in fact, the expansion 



tr[log(l-H2e,B)] 



fc=i 



k 



{2eq)Hv[B''],i75) 



valid for sufficiently small values of eg ^. Eq. I|74|l 
implies that only even powers of eg can con- 
tribute. This fact together with the observation 

^For instance, one can take 2\eq\ smaller than the number 
min [m,,l/||_B||2e,=m,]- 
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that B is an anti-Hermitian operator (B^ = —B, 
because 7 • V is anti-Hcrmitian) allows us to write 

~^ (-l)"+i 



i.h.s. of (Cgi 



n=l 



2n 



(2e,)2»tr[(BtB)»] 



1 



= -tr[log(l + Ae^B^B)] . 



(76) 



The last equality proves that tr[log(l + 2eqB)] is 
a non- negative quantity, implying that det[l + 
2eqB] is real and strictly positive. 

We conclude with two technical observations: 
1) since B is anti-Hermitian, its spectrum is 
purely imaginary, thus 1 + 2eqB cannot have 
any vanishing eigenvalue; 2) eq. (|76l) . which was 
proved for sufficiently small e,, can be extended to 
the actual physical value of the mass splitting by 
analyticity (see the chain of equalities in eq. H72|l 'l. 
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